CURVATURE BASED TRIANGULATION OF METRIC 
MEASURE SPACES 



EMIL SAUCAN 

Abstract. We prove that a Ricci curvature based method of triangula- 
tion of compact Riemannian manifolds, due to Grove and Petersen, ex- 
tends to the context of weighted Riemannian manifolds and more general 
metric measure spaces. In both cases the role of the lower bound on Ricci 
curvature is replaced by the curvature-dimension condition CD{K, N). 
We show also that for weighted Riemannian manifolds the triangula- 
tion can be improved to become a thick one and that, in consequence, 
such manifolds admit weight-sensitive quasimeromorphic mappings. An 
application of this last result to information manifolds is considered. 

Further more, we extend to weak CD(_ft', A'') spaces the results of 
Kanai regarding the discretization of manifolds, and show that the vol- 
ume growth of such a space is the same as that of any of its discretiza- 
tions. 



1. Introduction 

The existence of triangulation on geometric spaces can hardly be under- 
estimated both in Pm'e and in Apphed Mathematics, in particular that of 
certain special types (see, e.g. [5],[n],[12],[3I],[3Z],[3S],[Z],and[M],[2], 
[50] , [13], respectively). 

For Riemannian manifolds, a number of possible constructions exist, in- 
cluding those that produce special types of triangulations. Amongst them 
we mention, in chronological order and without any pretention of being ex- 
haustive, [H], PO], [31], [H] (Theorem 10.3.1), [3Zj, [3- 

We concentrate here on the method employed in [20]. (See also |3j for a 
history of this approach). The advantage of this method is, besides its ele- 
gance, the fact that it is highly geometric in spirit, using solely the intrinsic 
geometric differential properties of the manifold. The purpose of the present 
note is to show that the construction devised in [20] can be extended, almost 
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without any modifications to metric measure spaces (except, of course, the 
obvious necessary adaptations to the more general context). 

The reason behind this is dual: On one hand there exists a feeling in the 
community, that, although the tools and results developed by Gromov, Lott, 
Villani, Sturm and others are most elegant, they lack, so far, any concrete 
and efficient application. We wish, therefore, to further emphasize that the 
notions of curvatures for metric measure spaces are highly natural by giving 
an extension of a classical problem and its (also classical) solution, to the 
new context. On the other hand, it appears that there exists a real interest 
in the triangulation and representation of the information manifold, that is 
the space of parameterized probability measures (or the statistical model) 
equipped with the Riemannian metric induced by the Fischer information 
metric onto the Euclidean sphere (see Section 4.1.3 below). 

The reminder of this paper is structured as follows: As already mentioned 
above, we present, in Section 2, the construction of Grove and Petersen. 
(This and the following section are the most "didactic", thus we proceed 
rather slowly; however, afterwards the pace of the exposition will become 
more brisk.) Next, in Section 3, we introduce the necessary notions and 
results regarding curvatures of metric measure spaces. (Unfortunately, the 
style is rather technical, since we had to cover quite a large number of 
definitions and results.) Section 4 constitutes the heart of our paper, in the 
sense that we show therein how to extend Grove-Petersen construction to 
the metric measure setting. We follow, in Sections 5, with an extension, 
to weak CT){K, N) spaces, of Kanai's results regarding the discretization of 
manifolds, and show, in particular, that the volume growth of such a space 
coincides with that of any of its discretizations. In a sense, this represents 
the second part of the present paper, related to, yet distinct, from the main 
triangulation problem considered in the previous sections. We conclude, in 
Section 6, with a few very brief remarks. 

2. The Grove-Petersen construction 

In the following M"' = (M", g) denotes a closed, connected n-dimensional 
Riemannian manifold such that it has sectional curvature /cm bounded from 
below by k, diamM" bounded from above by D, and VolM"' is bounded 
from below by v. Since only volumes of balls arguments are employed, one 
can replace the last condition by the more general one RicM ^ (n — l)k (see, 
e.g., [32]). In fact, this very relaxation of the conditions actually helps us 
formulate the more general problem we are dealing with in the sequel. 

The basic idea is to use so called efficient packings: 

Definition 2.1. Let pi, . . . ,pnQ be points € M", satisfying the following 
conditions: 

(1) The set {pi, ■ ■ ■ ,Pno} is an e-net on M", i.e. the balls /3"(pfc,£)) 
/c = 1, . . . , no cover M"; 
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(2) The balls (in the intrinsic metric of M") (3^{pk,£/2) are pairwise 
disjoint. 

Then the set {pi, . . . ,Pno} is called a minimal e-net and the packing with 
the balls /3"'{pk,£/2), k = 1, . . . ,no, is called an efficient packing. The set 
{{k, I) \ k,l = 1, . . . ,nQ and /3'^{pk,£)r\f3^{pi,£) 7^ 0} is called the intersection 
pattern of the minimal e-net (of the efficient packing). 

Efficient packings have the following important properties, which we list 
below (for proofs see pO]): 

Lemma 2.2 (|20j. Lemma 3.2). There exists ui = ni{n,k, D), such that if 
{pi, . . . ,Pno} is a minimal £-net on M", then uq < ni. 

Lemma 2.3 (t2Qj, Lemma 3.3). There exists n2 = n2{n, k, D), such that for 
any x G M", |{j | j = 1, . . . , no and /3"(a;, e) fl j3'^{pj,£) 7^ 0}| < n2, for any 
minimal £-net {pi, . . . ,Pno}- 

Lemma 2.4 ([^Oj, Lemma 3.4). Let M", M2, be manifolds having the same 
bounds k = ki = k2 and D = Di = D2 (see above) and let {pi, . . . ,Pno} and 
{qi, . . . ,qno} be minimal £-nets with the same intersection pattern, on , 
M2 , respectively. Then there exists a constant = n^(n,k,D,C), such 
that if d{pi,pj) < C ■ £, then d{qi,qj) < ■ £. 

This properties suffice to provide us with a simple and efficient, even 
if crude, triangulation method of closed, connected Riemannian manifolds. 
Indeed, by using the properties above, one can construct a simplicial complex 
having as vertices the centers of the balls (3^{pk,£)., as follows: Edges are 
connecting the centers of adjacent balls; further edges being added to ensure 
the cell complex obtained is triangulated to obtain a simplicial complex. 
One can ensure that the triangulation will be convex and that its simplices 
are convex, by choosing e = ConvRad(M'^), where the convexity radius 
ConvRad(M"') is defined as follows: 

Definition 2.5. Let M"' be a Riemannian manifold. The convexity radius 
of M" is defined as inf{r > | f3'^{x, r) is convex, for all x G M"}. 

This follows from the fact that /?" (x, ConvRad(M")) C /S*" (x, InjRad(M")) , 
(since ConvRad(M") > ilnjRad(M") - see, e.g. [3]). Here InjRad(M") de- 
notes the injectivity radius: 

Definition 2.6. Let M" be a Riemannian manifold. The injectivity radius 
of Af" is defined as: InjRad(M") = inf{Inj(x) | x £ M"^}, where Inj(2;) = 
sup {r I e^Px\B'^ix,r) is a diffeomorphism}. 

Note that by a classical result of Cheeger (see, e.g. [3]), there is a universal 
positive lower bound for InjRad(M) in terms of k, D and v, where v is the 
lower bound for the volume of M. It is precisely this result (and similar 
ones - see also the discussion below) that make the triangulation exposed 
above a simple and practical one, at least in many cases. 
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Remark 2.7. Lemmas 2.2 - 2.4 above represent part of the tool£| employed 
in the proof of main theorem of [20], namely: 



Theorem 2.8 ([20j, Theorem 4.1). Let M^,M2 he two manifolds having the 
same upper diameter hound D, as well as the same lower hounds k and 
V, on their curvatures and volumes, respectively. Then there exists e = 
e{n,k, D,v) such that, if Mi and M2 have minimal packings with identical 
intersection patterns, then they are homotopy equivalent. 

Unfortunately, the condition regarding the sectional curvature bound can- 
not be replaced easily by a similar one regarding Ricci curvature, not even 
in the classical Riemannian casell and evidently not in the more general 
setting adopted in this paper. 

We bring, for reference, the proofs of the Lemmas above, and we do this 
almost verbatim: 

Proof of Lemma 2.2. Let {pi, . . . ,Pno} be a minimal e-net on M" and let p 
be a point in - the fe-space form. Then, by the classical Bishop-Gromov 
Theorenil 

\o\B{p,r) \o\B{p,r) 

\o\B{p, R) - \oW{p, R)' ' 

for any p € Af". 

Let io such that Voli?(pjQ, e/2) is minimal. By 2.1.(2) it follows that 

VolM" \o\B{p,D) 
no<TTT7T7 ^< 



\oW{pi„e/2) - \o\B{p,e/2) ' 
(To obtain the last inequality, just take, in Bishop-Gromov Theorem, R 
diamiif" <D.) 

The desired conclusion now follows by taking 



notable part of the other tools being represented by a generalization of Cheeger's 
"butterfly" construction [lUj (see also I15j). 

See Berger |3] for a brief discussion on the results employing Ricci curvature bounds. 
■^It basically states that the volume of balls in a complete Riemannian manifold {M" , g) 
satisfying Ric > (n — l)k, does not increase faster than the volume of balls in the model 
space form, more precisely that, for any x £ M = M" , the function 

VolB(x,r) 



is nonincreasing (as function of r), where 



rin 
OK 



K(r) 



UK >0 
ii K = 
sinh,/^r) ifi4:<0 
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YolB{p,D) 
YolB{p, e/2) 



□ 



Proof of Lemma 2.3. Let ji, js he such that B{x,e)r\B{pj-,e) / 0. Then 
B{pj^,e/2) cB{x,5e/2). 

Let k € {1, ...,s} be such that B{pj^^£ /2) has minimal volume. Then (as 
in the proof of Lemma 2.2) it follows that: 

^ Vol^(x,5g/2) ^ Volg(pj^,9g/2) ^ Vol^(p, Qe/2) 
^- YolB{pk,e/2) - Vols (p,,, e/2) " VolS(p,e/2) ' 
where p is as in the proof of the previous lemma. But 

YolBip,9e/2) _ f^'^' S^K{r)dr 



\o\B{p,e/2) jl/^Sl{r)dr ' 

and the function 

lf'Sl{r)dr 

extends to a continuous function h : [0,1?] — > M+, (since h{e) — )• when 
e ^ 0). 

□ 

Remark 2.9. It is important to note that 712 is independent of e. 

Proof of Lemma 2.4- Evidently, since d{pi,pj) < C ■ e, it follows that pj G 
B{pi,C -e). Thus, precisely as in the proof of the previous lemma, it follows 
that there exists n' = n'{C), 

n [C) = max — -j- , 

rJ"Sl{r)dr 

such that at most n' of the balls B{pi,e/2), B{pn' ,£/2) are included in 
Bip,,{C + ^)e/2. 

Since {pi, . . . ,Pno} and {qi, . . . , qng} have the same intersection pattern, 
it follows that d{qi,qj) < n3(C), where n3(C) = 2[n'(C) — 1]. 

□ 

3. RiCCI CURVATURE OF METRIC MEASURE SPACES 

We bring here the definitions and results necessary in the following sec- 
tion. However, since proofs are not needed, we omit them and send to the 
source (usually [33] )j^ 

We begin with the basic (and motivating) generalization, namely: 



^However, the reader should be aware that, while we mostly refer to [4j for convenience, 
many (if not all) the results therein appeared first in [27] and/or |41| . 
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3.1. Riemannian manifolds. Let M = M" be a complete, connected n- 
dimensional Riemannian manifold. 

One wishes to extend results regarding Ricci curvature to the case when 
M" is equipped with a measure that is not dVol. Usually (at least in our 
context) such a measure is taken to be of the form 

u{dx) = e-^(^')Vol(dx) , 
where V : M" — t- M, F G C^(M). Note also that any smooth positive 
probability measure can be written in this manner. Then {M,d,i'), where 
d is the geodesic distance, is a metric measure space. 

Remark 3.1. A standard measure in the context of Image Processing (but 
not only) would be the gaussian measure on M": 

(n) ^ e-I^I'drc 
^ (27r)"/2 ■ 

However, more realistic measures can (and, indeed, should) be considered 
for imaging purposes - see [21] , [23] . 



To preserve geometric significance of the Ricci tensor, one has to modify 
its definition as follows: 

(3.1 RicAT^ = Ric + VV 

N — n 

Here (S" W is a quadratic form on TM'"" , and V'^V is the Hessian 
matrix Hess, defined as: 

{vv vv)^{v) = (vy(x) • vf . 

Therefore 

Ric;v,.(7) = (Ric + W)(7) - . 
Here is the so called effective dimension and is to be inputed. 

Remark 3.2. (1) If A < n then RicTv.i/ = —oo 

(2) If A = n then, by convention, x oo = 0, therefore ()3.ip is still 
defined even if W = 0, in particular Ric„ voi = Ric (since, in this 
case y = 0). 

(3) If A = oo then Ricoo,,. = Ric + W. 

The Ricci curvature boundedness condition of the classical Bishop-Gromov 
is paralleled in the case of smooth metric measure spaces by the following 
(rather obvious) 

Definition 3.3. (M, d, v) satisfies the curvature- dimension estimate CD(i^, A^) 
iff there exist X € M and A^ € [1, oo], such that RicTv^j^ > K and n < N . (If 
V = dVol, then the first condition reduces to the classical Ric > K.) 
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Remark 3.4. Intuitively, "M has dimension n but pretends to have dimen- 
sion N. (Identity theft)"| 

The need for such a parametric dimension stems, in particular, from the 
desire to extend the Bishop-Gromov Theorem to metric spaces (or more 
precisely, to length spaces), for which no innate notion of dimension exists. 

Remark 3.5. For a number of equivalent conditions, see [44], Theorem 14.8. 

Theorem 3.6 (Generalized Bishop-Gromov Inequality). Let M be a Rie- 
mannian manifold equipped with a reference measure u = e~^Vol and sat- 
isfying a curvature- dimension condition CD(K, N), K < N < oo. 
Then, for any x € M , the function 

u[B{x,r)] 
""^'^ j;S^{t)dt' 
is nonincreasing ( as function of r ), where 



N-l 



sin^/j^t] ifK>0 



t^'' ifK = 



N~l 



(sinh^j^tj ifK<0 



Proof See [13], p. 499-500. □ 

3.2. Metric Measure Spaces. To further extend the Bishop-Gromov The- 
orem, we have first to introduce a proper generalization of "Ricci curvature 
bounded from below" . For this, quite a number of preparatory definitions 
are needed, that we briefly review here. 

Definition 3.7. Let {X,fi) and (Y,!^) be two measure spaces. A coupling 
(or transference (transport) plan) of fi and is a measure vr on X x y 
with marginals fi and (on X and Y, respectively), i.e. such that, for all 
measurable sets A C X and B C Y, the following hold: 7r[A x y] = fj,[A] 
and tt[X x B] = u[B]. 

Definition 3.8. Let {X, n) and {Y, u) be as above and let c = c(x, y) be a 
(positive) cost function on X x y. Consider the Monge-Kantorovich mini- 
mization problem: 



(3.2) inf / c{x,y)dTT{x,y) , 

JxxY 

where the infimum is taken over all the transport plans. The transport plans 
attaining the infimum are called optimal transport (transference) plans. 



5 J. Lott [26]. 
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Before we can proceed, we must recall the following definition and facts: 

Definition 3.9. Let {X, d) be a Polish space, and let P{X) denote the set of 
Borel probability measures on X. Then the Wasserstein distance (of order 
2) on P{X) is defined as 

(3.3) W2i^^, u) = (^inf ^ d((x, y)^d7r{x, y) 

where the infimum is taken over all the transference plans between fj. and v. 
Definition 3.10. The Wasserstein space P2{X) is defined as 

(3.4) P2{X) = G P{X) I j dixo,x)^fi{dx) < oo} , 

where Xq (z X is an arbitrary point. 

The definition above does not depend upon the choice of xq and W2 is a 
metric on P2{X). 

Facts 3.11. Let {X,d) and P2{X) be as in the definition above. Then 

(1) P2iX) is a Polish space. 

(2) If X is compact, then P2iX) is also compact. 

Definition 3.12. Let {X,d) be a compact, geodesic, Polish space and let 
r = {7 : [0,1] — > X I 7 a minimal geodesic}, and denote by et : F ^ X 
the (continuous) evaluation map, et : (7) = 7(i). Let E : T ^ X x X he 
defined as £'(7) = (69(7), 61(7)). A dynamical transference plan is a pair 
(7r,n), where vr is a transference plan and n is a Borel measure, such that 
ii^^n = vr. (vr, H) is called optimal if vr is optimal. 

Definition 3.13. Let H be an optimal dynamical transference plan. Then 
the one-parameter family {/^t}te[o,i] , = (et)#n is called a displacement 
interpolation 

We can now quote the following result ([44], Theorem 7.21 and Corollary 
7.22), connecting the geometry of the Wasserstein space to classical mass 
transport: 

Proposition 3.14. Any displacement interpolation is a Wasserstein geo- 
desic, and conversely, any Wasserstein geodesic is obtained as a displace- 
ment interpolation from an optimal displacement interpolation. 

Definition 3.15. Given N G [l,oo], the displacement convexity class T>Cm 
is defined as the set of convex, continuous functions U : M+ — )• M, ?7 G 
C2(M+ \ {0}), such that C/(0) = and such that 

rU'{r) - U{r) 
is nondecreasing (as a function of r). 
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Remark 3.16. For equivalent defining conditions for the class "DCn see [Hj . 
Definition 17.1. 

Definition 3.17. Let (X, d, i/) be a a locally compact metric measure space, 
such that the measure v is locally finite, and let U he a continuous, convex 
function U : R+ ^ R, U e C^(R+ \ {0}), such that U{0) = 0. Consider 
a measure fj, on X, having compact support, and let /j, = pu + fig he its 
Lebesgue decomposition into absolutely continuous and singular parts. 

Then we define the (integral) functional Ui^ (with nonlinearity U and 
reference measure v) by 



Moreover, if {'i:{dy\x)}x&x is a family of probability measures on X and 
if /3 : U X U ^ (0, oo] is a measurable function, we define an (integral) 
functional U^^,^ (with nonlinearity U, reference measure u, coupling it and 
distortion coefficient (3) by: 



Usually (e.g. in the definition of weak CD{K^ N) spaces) /3 is taken to be 
the reference distortion coefficients: 

Definition 3.18. Let x,y be two points in a metric space (X, d), and con- 
sider the numbers K €, N £ [l,oo] and t € [0, 1]. We define the reference 
distortion coefficients /3|^'^'*(x, y) as follows: 



(3.5) 




(3.6) 




(1) If t G (0, 1] and 1 < iV < oo, then 



+ 00 



if K > and a > vr 




if K > and a € [0, vr] 



(3.7) 



{x,y) 



= < 



1 



if K = 




where 



(3.8) 




(2) In the limit cases — )• 1 and — )• oo, define 
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(3.9) 

and 



+00 if > , 
1 if if < ; 



(3.10) (3f'^\x,y) = eT(^-''^'^(-'y)- 
(3) If i = 0, then 

(3.11) /3f'^)(x,y) = l. 

Remark 3.19. If X is the model space for CD{K,N) (see [S] p. 387), then 
pi^'^) ig the distortion coefficient on X. 

We can now bring the definition we are interested in: 

Definition 3.20. Let {X,d,i') be a locahy compact, complete, fj-finite 
metric measure geodesic space, and let if G M, G [l,oo]. We say that 
{X, d, v) satisfies a weak CD(if , A'^) condition (or that it is a weak CD(iC , A^) 
space) iff for any two probability measures /^Oi^i with compact supports 
Supp/ii, Supp/i2 C Suppi^, there exist a displacement interpolation /^to<t<i 
and an associated optimal coupling vr of fJ.Q,fJ-i such that, for all U G VCn, 
and for all t € [0, 1], the following holds: 

(3.12) U,{^lt) <{l-t) ifio) + 1 c/f ^ (/^i) 
(Here we denote vr = S^n, where S{x,y) = {y,x).) 

Remark 3.21. In fact, the geodesicity condition is somewhat superfluous, 
since a locally compact, complete metric space is geodesic (see, e.g. [2T] , 
9.14). 

Fortunately, a version of the Bishop-Gromov Inequality also holds in gen- 
eral metric measure spaces; more precisely we have 

Theorem 3.22 (Bishop-Gromov Inequality for Metric Measure Spaces). 
Let [X^d^v) he a weak CD{K,N) space, N < 00, and let xq G Suppz/. 
Then, for any r > 0, u[B{xQ),r] = u[B[xQ\,r]. Moreover, 

(3.13) "[''[""l''] 



is a nonincreasing function of r. 

Proof See LMJ, p. 806. □ 
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For US, the following corollary of this extension of the Bishop-Gromov 
Theorem is very important, since it allows for a lower hound on the number 
of balls in a efficient packing: 

Corollary 3.23 (Measure of small balls in weak CD{K, N) spaces). Let 

{X,d,iy) be a weak CD{K,N) space, N £ [l,oo), and let z G Suppz^. Then, 
for any R > 0, there exists c = c{K,N,R) such that, if B{x,r) C B(z,R), 
then: 

(3.14) iy[B{x, r)] > {cu[B{z, R)])r^ . 

Since, by [H], Theorem 29.9, any smooth CT){K, N) metric measure space 
is also a weak CT){K, N) space, it follows that a lower bound on the number 
of balls in an efficient packing also exists for smooth metric measure spaces. 

Moreover, an expected doubling property also holds: 

Corollary 3.24 (Weak CD{K, N) spaces are locally doubling). Let {X, d, v) 

he a weak CD{K,N) space, N G [l,oo), and let z G Suppz^. Then, for 
any fixed ball B{z,R) C X, there exists C = C{K,N,R) such that, for all 
r G {0,R), i'[B{z,2r)] < Ci'[B{z,r)]. In particular, if diamX < D, then X 
is C -doubling, where C = C{K,N,D). 

4. Triangulating metric measure spaces 
4.1. Smooth metric measure spaces. 

4.1.1. The basic construction. The lemmas regarding efficient packings trans- 
late to the context of smooth metric measure spaces with little, if any, mod- 
ifications: 

Lemma 4.1 (Lemma 2.2 on manifolds with density). Let {M^ ,d,u), v = 
e~^Vol, be a compact(closed) smooth metric measure space such that Ric^^y > 
K, for some K < N < oo, and such that diamM" < D. Then there 

exists ni = ni{N, K, D), such that if {pi, . . . ,Pno} is a minimal e-net on 
M", then hq <ni. 

Proof. Since the only essential ingredient of the proof of Lemma 2.2 is the 
Bishop-Gromov volume comparison theorem, and since by Theorem 13.61 its 
analogue also holds for smooth metric measure spaces, the proof follows 
immediately precisely on the same lines as that of Lemma 2.4. □ 

Lemma 4.2 (Lemma 2.3 on manifolds with density). Let {M^,d,u) he as 
in Lemma \4.1\ There exists n2 = n2{N,K,D), such that, for any x G , 
\{j |j = 1, . . . , no and /3"(x, e) PI /3^{pj,e) ^ 0}| < n2, for any minimal e-net 
{Pi, ■ ■ ■ ,Pno}- 

Proof. Again, as in the proof of the previous lemma, the only relatively new 
ingredient is the generalized Bishop-Gromov Theorem, which as we have 
seen, holds for smooth metric measure spaces. □ 
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Lemma 4.3 (Lemma 2.4 on manifolds with density). Let (M^ ,di,i^i) and 
{M2,d2,i'2) be as in Lemma\4^ and let {pi, . . . jPn^} and {qi, . . . , Quq} be 
minimal e-nets with the same intersection pattern, on M", M2 , respectively. 
Then there exists a constant = n3{N,K,D,C), such that if di{pi,pj) < 
C ■ e, then d2{qi,qj) < ■ e. 

Proof. Idem, the proof follows along the lines of the original proof (of Lemma 
2.4), by force of Theorem 13.61 □ 

Remark 4.4. Note that the constants ni,n2,n3, as functions of the effective 
dimension N, rather than the topological one n, are only very weekly de- 
pendent (in general) on the geometry of the manifold M" , via the inequality 
condition n < N . While they still will guarantee the existence of minimal 
e-nets with the required properties, the metric density of this nets will be 
different from the purely geometric one given by the classical Grove-Petersen 
construction, hence so will be the shape (or "thickness", see Definition 4.8 
below) of the simplices of the resulting triangulation. 

The basic triangulation process now carries to the context of smooth 
metric measure spaces without any modification. The convexity of the tri- 
angulation follows exactly like in the classical case, since the injectivity and 
convexity radius are solely functions of the metric d and not of the weighted 
volume V. 

Remark 4.5. As we have already noted in Remark 2.7, the lower bound on 
the sectional curvature is essential for the proof of the Homotopy Theorem 
12.81 Therefore, one cannot formulate a similar theorem for the weighted 
manifolds, without imposing the additional constraint on sectional curva- 
ture. The original argument of Grove and Petersen still holds for the Rie- 
mannian structure, while measure has no other role than determining the 
existence and metric density of the triangulation vertices. In consequence, 
we can only obtain a rather weak result, that we bring here (for complete- 
ness' sake) only as corollary, namely: 

Corollary 4.6 (Theorem l2.8l for smooth metric measure spaces). Let (MJ^, di^vi), 
{M2 ,d2,i'2) , Vi = e~^'dVol, Vi G C^(M), i = 1,2 be smooth, compact metric 
measure spaces satisfying CD{K, N) for some X G R and \ < N < 00, 
and such that diamM" < D, VolM" < v, i = 1,2 and, moreover, having 
the same lower bound k on their sectional curvatures. Then there exists 
e = £{N, K, k, D, v) such that, if Mf, M2 have minimal packings with iden- 
tical intersection patterns, they are homotopy equivalent. 

4.1.2. Thick triangulation and quasimeromorphic mappings. We can strengthen 
the simple result above, to render a "geometrically nice" triangulation, 
namely we can formulate the following 

Proposition 4.7. Any smooth, compact metric measure space {M'^,d,v) 
satisfying CT)[K, N) admits a (p* -thick triangulation, where (p* = 'p*{n, d, v). 
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Recall that thick (or fat) triangulations are defined as follows: 

Definition 4.8. Let r C M" ; < A; < n be a fc-dimensional simplex. The 

thickness (/? of r is defined as being: 



The infimum is taken over all the faces of r, a < r, and Yolj{a) and diamo" 
stand for the Euclidian j'-volume and the diameter of a respectively. (If 
dimo" = 0, then Volj(c7) = 1, by convention.) A simplex r is ipo-thick, 
for some (po > 0, if (/'(''") ^ ^^o- A triangulation (of a submanifold of M") 
T = {(Ti}i£i is ipo-thick if all its simplices are (/jg-thick. A triangulation 
T = {cTj}igi is thick if there exists > such that all its simplices are 
ipQ-thick. 

(The definition above is the one introduced in [llj. For some equivalent 
definitions of thickness, see [8], [28], [3l], [l2], [l3].) 

Remark 4.9. By [11], pp. 411-412, a triangulation is thick iff the dihedral an- 
gles (in any dimension) of all the simplices of the triangulation are bounded 
away from zero. 

Proof. Since the geometry of the manifold is not affected by the presence of 
the measure u, the 6-transversality and e-moves arguments of Cheeger et al. 
apply unchanged - see [11] for the lengthy technical details. 

Moreover, the resulting triangulation is "z^-sensitive" , insomuch as the 
metric density of the vertices, hence the thickness if* of the triangulation's 
simplices is a function of the measure v (as shown by Lemmas l4.1H4.3l above). 



Remark 4.10. As the proof above shows, the role of the measure in deter- 
mining the thickness of the triangulation is somewhat marginal. Therefore, 
it would be useful to modify the thickness condition so that it will reflect 
(and adapt to) the presence of the measure v. (This may prove to be even 
more relevant in the case of weak CD(X, A^) spaces.) However, if the new 
angles /-k,n depend continuously on the standard ones, then, given that 
V = e~^Vol and V G C^(M), the continuity and compactness arguments 
involved in the proof of the classical case also hold for generalized (dihedral) 
angles, a generalization of Proposition 14.71 will follow immediately (as will, 
in consequence, the corresponding one of Corollary 14.131 below) . 

Remark 4.11. For Riemannian manifolds, on can construct thick triangula- 
tions, for non-compact manifolds witlH or without boundary, (see, e.g. [31] . 
[37| . [39]). However, for weighted manifolds this is far from obvious, because 
there exists no a priori information regarding the properties of the restric- 
tion of the measure v to the (n — l)-dimensional manifolds required in the 
process of the triangulation. 



(4.1) 



ip = ip{T) = inf 



dim cr = j 




□ 



at least for many cases 
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Before proceeding further, we have to remind the reader the definition of 

quasimeromorphic mappings: 

Definition 4.12. Let M'^,N'^ be oriented, Riemannian n-manifolds. 

(1) / : M" N"^ is called quasiregular (qr) iff 

(a) / is locally Lipschitz (and thus differentiable a.e.); 
and 

(b) < |/'(x)|" < KJfix), for any x G M"; 

where f'{x) denotes the formal derivative of / at x, |/'(x)| = sup \f'{x)h\, 

\h\ = 1 

and where J/(x) = detf'{x); 

(2) quasimeromorphic (qm) iff = S", where is usually identified 
with = U {oo} endowed with the spherical metric. 

The smallest number K that satisfies condition (b) above is called the outer 
dilatation of /. 

Using the arguments of, say, [37] we obtain the following 

Corollary 4.13. Any smooth, compact metric measure space (M^jdjiy) sat- 
isfying CD{K,N) admits a non-constant quasimeromorphic mapping f : 

4.1.3. An Application - The Information Manifold. The method of triangu- 
lation and quasimeromorphic mapping of weighted Riemannian manifolds 
presented above represents a generalization of a result classical in informa- 
tion geometry. To be somewhat more concrete we present it briefly here (for 
more details see, e.g. [l], [T7]). 

Let A be a finite set, let fi{x), i = 1,2 be bounded distributions on A, 
and let Pi(x) = v^^^^^tttt; viewed as probability densities on A. The relative 
information between pi and p2 (or the Kullback-Leibler divergence) is defined 
as 

KL{pi\\p2) = ^pilog 

A 

represents a generally accepted measure of the divergence between the two 
given probabilities, but, unfortunately, it fails to be a metric. However, it 
induces a Riemannian metric on P{A) - the manifold of probability densities 
on A, namely the Fisher information metric: 




(4.2) 5Fischcr,p(A) = KL{p,p + A) = V -^- , 

where p G -P(^) is given and A represents an infinitesimal perturbation. 
(Here, to ensure that p + A will also be a probability density, the following 
normalization is applied: A{x)dx = 0.) 
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Remark 4.14. It turns out (see [T]) that the Fisher information can be 
written as Riemannian metric (in standard form) in the following form: 
SFischer,- = (gij), where 



where 9^,... jd'^jk = \A\ represent the coordinates on ctq, 1 = logp is the 
so called log-likelihood and Ep(fg) denotes the expectation of fg, Ep(fg) = 

J fgdp- 

The correspondence p{x) ^ u{x) = 2y^p(x) maps the probability sim- 
plex (To = {p{x) \ x € A,p{x) > 0,^^p(x) = 1}, onto the first orthant of 
the sphere S = YIa '^{x)'^ = 4. This mapping preserves the geometry, in the 
sense that the geodesic distance between p,q G ctq), measured in the Fisher 
metric, equal the spherical distance between their images (under the map- 
ping above). Moreover, geodesies are mapped to great circles. (For details 
and more geometric insight regarding the probability simplex and the map 
above, see [l].) 

To summarize: The quasimeromorphic mapping of a weighted Riemann- 
ian manifold {M^,d,v) onto the n-dimensional unit sphere S", represents a 
generalization of the considerents above in two manners: 

(a) It allows for the mapping with controlled and bounded distortion 
(i.e. qm) of a more general class of Riemannian manifolds (with arbitrary 
metrics) endowed with a variety of (probability) measures, and not just of 
the standard statistical model; 

(b) It permits the reduction to the study of the (geometry of the ) stan- 
dard simplex in of the geometry of the whole information manifold, and 
not just of the probability simplex. (As far as application of such gener- 
alizations are concerned we mention here only those related to signal and 
image processing (see e.g. [16], [IT], respectively), and information theory 
(see, e.g. [I]).) 

To be sure, a number of minor technical details need to be considered, 
namely the fact that we allow for multidimensional spheres, i.e. we consider 
also multivariate distributions, as opposed to the simpler notation adopted 
above; and the use of the standard simplex (first orthant) in the unit sphere, 
rather than in the sphere of radius 4, which is trivial, giving that dilation is 
a (elementary) conformal mapping. 

Remark 4.15. As we have seen above, bounds on the injectivity (hence con- 
vexity) radius exist in terms of curvature (in conjunction with diameter and 
volume), due to what represent, by now, classical results. This conducts us 
to formulate the following question: 

Question 1. Given a smooth metric measure space, is it possible to find 
hounds on the injectivity radius in terms of the CD{K,N) condition? 



(4.3) 
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In fact, it should be even possible to find bounds on the conjugate ra- 
dius (see, e.g. [3j), since conjugate points appear when the volume density 
vanishes (see |32j . Section 2.4.4) and this is a function of the determinant 
Jacobian, which is controlled by the Ricci curvature (see [32], Sections 2.4.2 
and 6.3.1, fii]). 

4.2. Weak CD{K, N) spaces. Using Theorem 13. 221 the proofs of the basic 
lemmas, for metric measure spaces, follow immediately. We do, however, 
bring below their statements, for the sake of completeness: 

Lemma 4.16 (Lemma 2.2 on weak CT)(K,N) spaces). Let (X, d, z/) be a 

compact weak CD{K, N) space, N < oo, such that Suppj^ = X and such that 
diamX < D. Then there exists ui = ni(K, N, D), such that if {pi, . . . ,Pno} 
is a minimal e-net in X, then uq <ni. 

Remark 4.17. Note that, since < oo, the condition Suppi^ = X imposes 
no real restriction on X (see [53], Theorem 30.2 and Remark 30.3). 

Lemma 4.18 (Lemma 2.3 on weak CD{K,N) spaces). Let {X,d,u) be a 
compact weak CD{K, N) space, N < oo, such that Suppz^ = X and such that 
diamX < D. Then there exists n2 = n2{N,K,D), such that, for any x G 
, \{j |j = 1, . . . , no and /3"(x, e) PI /3^{pj,e) ^ 0}| < n2, for any minimal 
e-net {pi, . . . ,Pno}- 

Lemma 4.19 (Lemma 2.4 on weak CD{K, N) spaces). Let (M", di, ui) and 
(M^, d2, 1^2) be as in Lemma [4.16\ and let {pi, . . . ,Pno} o,nd {qi, . . . , he 
minimal e-nets with the same intersection pattern, on M", M2 , respectively. 
Then there exists a constant = n3{N, K, D,C), such that if di{pi,pj) < 
C ■ e, then d2{qi, Qj) < ■ e. 

The existence of the triangulation follows now immediately, since, by 
definition, weak CD{K,N) spaces are geodesic. Moreover, in nonbranching 
spaces, the geodesies connection two vertices of the triangulation are unique 
a.e. (see j44j . Theorem 30.17). Recall that a geodesic metric space X is 
called nonbranching iff any two geodesies 71, 72 : [0, t] ^ X that coincide on 
a subinterval [0,to]iO <tQ <t, coincide on [0,t]. 

However, it is hard to ensure the convexity of the triangulation. Sadly, 
many weak CT){K, N) spaces of interest fail to be locally convex. Fortu- 
nately, local convexity does hold for an important class of metric measure 
spaces: Indeed, by [34], Alex[ir] C CD((m — l)K,m), where Alex[if] de- 
notes the class of m-dimensional Alexandrov spaces with curvature > K, 
equipped with the volume measured But, by [33,], Lemma 4.3, any point in 
Alexandrov space of curvature > K has a compact neighbourhood which is 
also an Alexandrov space of curvature > K. 

Remark 4.20. The result of Perelman and Petrunin is, in fact, a bit stronger, 
in the sense that the diameter of the said compact neighbourhood is speci- 
fied. On the other hand, their proof uses the rather involved (even if by now 



In fact, Petrunin provides the full details of the proof only for the case K — 0. 
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standard) tool of Gromov-Hausdorff convergence (of Alexandrov spaces), 
amongst others. 

An alternative, more intuitive (at least in dimension 2) proof can be given, 
using the notion of Wald-Berestovskii curvature (see [4J, ^SSj). However, 
since we do not want to further encumber the reader with more definitions, 
and since Alexandrov spaces represent by now a widely accepted tool in 
Geometry, that generates important research on its own, we shall bring the 
full details of this proof, and of the whole construction, elsewhere. 

Remark 4.21. Since the balls f3^{pk,£) cover M^, any sequence V{£m)-, of 
efficient packings such that — >• when m — )• oo, generates a discretization 
{Xm, d^Vrn), in the sense of |6]: Consider the Dirichlet (Voronoi) cell complex 
(tesselation) = {C(pm,fc)fc} of centers pfc^m. and atomic masses Vm{Pm,k) = 
^[C{pm,k)\- Then taking Xm = {Pm,k\-, d the original metric of X and Um 
as defined, provides us with the said discretization. 

It follows, by [6], Theorem 4.1, if Vol(M") < oo, the sequence {Xmi d, Um) 
converges in the W2 (see [2], [H]) to a metric measure space and, moreover, 
if {X,d,i') is a weak CD{K,N) space, then, for small enough e, so will be 
{Xm, d, Um), but only in a generalized ( ^^rough" ) sense. (For details regarding 
the precise definition of rough curvature bounds and the proof of this and 
other related results, see [6j.) 

5. Discretizations and volume growth rate 

The compactness condition imposed in the previous sections stems from 
the need for finding estimates in Lemmas 2.2-2.4 (and their respective gen- 
eralizations) and as a basic requirement for Theorem 2.8 to hold. However, 
the basic geometric method employed for obtaining efficient packings holds 
for noncompact manifolds, as well (given a lower bound on the injectivity 
radius). 

Therefore, we discard in this section the compactness restriction, and 
concentrate, following Kanai [25] (see also [9]) on those properties of e-nets 
that hold also on unbounded manifolds, and mainly on the volume growth 
rate. 

Amongst these properties, the most basic is the one contained in condition 
(2) of Definition 2.1, that is that if Pi,Pj & M,i ^ j, then d{pi,pj) < 2e. We 
also impose the additional condition that the set J\f, satisfying the condition 
above, be maximal with respect to inclusion. 

Following Kanai [25], we call the graph G{N') obtained as in Section 
2 above (i.e the 1-skeleton of the simplicial complex constructed therein) 
a discretization of X, with separation e (and covering radius e) (or a e- 
separated net). Further more, we say that G{J\f) has bounded geometry iff 
there exists po > 0, such that p{p) < pQ, for any vertex p G Af, where p{p) 
denotes the degree of p (i.e. the number of neighbours of p). 

Remark 5.1. Note that Kanai's discretization is different from the one con- 
sidered in p. Evidently, there exists a close connection between the two 
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approaches: Instead of the combinatorial length on the edges of the graph 
G{M), consider the length of the geodesies (between the adjacent vertices). 
If one restricts himself to smooth metric measure spaces, then this graph can 
actually be embedded in M" (see remark after Definition 2.5). The Voronoi 
cell complex construction follows, of course, without any change. (One can 
adopt even a semi-discrete approach: Using combinatorial lengths for the 
edges of the graph and atomic measures (weights) for the vertices (i.e the 
nodes of the graph) equal to the measure of the corresponding Voronoi cell.) 

Moreover, the graph G considered above, together with the geodesic met- 
ric dg induced by the metric of M", represents a geodesic metric space 
and, as such, it possess, at each vertex v, a (not necessarily unique) metric 
curvatur^ k{v) (e.g. the Wald-Berestovskii curvature mentioned above). 
Therefore, a natural question arises in connection with the discretization 
discussed in Remark 4.210 

Question 2. What is the relation (for e small enough) between above, 
and the rough curvature bound of the corresponding space {X^md^Vm) from 
Remark 4.21? 

Remark 5.2. Note that in ^25j a more general covering radius is consid- 
ered. For the cohesiveness of the paper we have used e as covering radius. 
However, the proofs of all results using covering radius hold for any R. 

We start (presenting the results) with the following lemma: 

Lemma 5.3. Let {X,d,v) be a (weak) CT){K,N) space, K < 0, N < oo, 
and let M be a e-separated net. Then 

(1) 



\MriB{x,r)\ < 



pip) < 



for any x G X and r > 0. 



Proof. Bv Lemma l4.18l and Corollarv [3^231 {MnBix. r)| is finite. The precise 
bounds (1) and (2) are obtained in the course of the proofs of l4.19l and I4.18| 
(or rather in their classical versions) respectively. □ 

Definition 5.4 (Rough isometry). Let {X,d) and {Y,6) be two metric 
spaces, and let f : X ^ Y (not necessarily continuous). / is called a 
rough isometry iff 



^See 0]. 

^ See also [6], Example 4.4 and Section 5. 



TRIANG-MMSP 19 

(1) There exist o > 1 and 6 > 0, such that 

-d{xi,X2) - b< d{f{xi), f{x2)) < ad{xi,X2) + b , 
a 

(2) there exists ei such that 

U B{f{x),ei>0)=Y; 

(that is / is ei-full.) 

Remark 5.5. (1) Rough isometry represents an equivalence relation. 
(2) If diam(X), diam(y) are finite, then X, Y are roughly isometric. 

The basic result of this section is the following generalization of |25j 
Lemma 2.5 (see also [9j, Theorem 4.9). 

Theorem 5.6. Let {X,d,i') be a weak CD{K,N) space and let G be a 
discretization of X. Then {X,d) and (G, d), where d is the combinatorial 
metric, are roughly isometric. 

The proof closely follows the one for the classical case, except for a nec- 
essary adaptation: 

Proof. By the construction of A/" and G{J\f), the following evidently holds: 

(5.1) d{pi,p2) < 2ed(pi,p2), Pi,P2 G A/", 

which gives the required lower bound (with a = l/2e and without b) without 
any curvature constraint. 

To prove the second inequality, one has, however, to assume that {X, d, u) 
is a weak CD{K,N). 

Then, given pi,p2 £ AA, chose a minimal geodesic 7 C G{J\f) with ends 
pi and p2. Put T.y = {y £ X \ d{y, 7) < e} and denote AA-y = M Ci T^. Then 
there exists a path in G{N) connecting pi and p2 and contained in AA-y (see 
[9], p. 196). Therefore 

d(jpi,P2) < \^f^\ . 

Let m be such that (m — l)e < d{pi,p2) < e, and let pi = qo, . . . , Qm = P2 
equidistant points on 7. Then 

,/ N d{pi,p2) 

d{qi-i,qi) = < £. 

m 

From Lemma 5.2 it follows that 

\My\r\B{qi,2£) <C = G{K,N,e)). 

Moreover, 
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m 
1=0 

We obtain, after a few easy manipulations that 
where Co,Ci are constants. 

□ 

Remark 5.7. Holopainen |22j has shown that if one only requires that d{pi,p2) < 
3e, instead of < 2e, as in Kanai's original construction, then the theorem 
above can be proved by purely metric methods, without any curvature con- 
straint. However, we have included here Kanai's version, for its Ricci cur- 
vature "flavor". Also, as already noted above, the construction is, in this 
case, geometrically natural. 

By Remark 5.3 (2) above, immediately follows 

Corollary 5.8. Any two discretizations of a weak CD{K,N) space are 
roughly isometric. 

Before proceeding further, we need to remind the reader the following 
definition: 

Definition 5.9 (Volume growth). For x € X and r > we denote the 
"volume" growth function by: 

(5.2) Vix,r) = i^[B{x,r)]. 

We say that X has exponential (volume) growth iff 

1^1 X r) 

(5.3) lim sup ^ ' ^ > , 

r— >-oo r 

and polynomial (volume) growth iff there exists k > such that 

(5.4) V{x,r) <C ■r'' , 
(C = const.) for all sufficiently large r. 

For the vertices of G{J\f), there exist two natural measures: the counting 
measure dr. (it(7W) = \^A\, for any A4 C J\f; and the "volume" measure dY: 
d'V{p) = p{p)di{p). However, in our context, one can substitute di for dV 
(and, indeed, for any measure /i absolutely continuous with respect to di) - 
see [9], page 197 - so we shall work with the counting measure, because its 
intuitive simplicity. 

Lemma 5.10 (Roughly isometric graphs have identical growth rate). Let G 
and T be connected, roughly isometric graphs with bounded geometry. Then 
G has polynomial (exponential) growth iff G has polynomial (exponential) 
growth. 
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Proof. See [25j, p. 399. □ 

Theorem 5.11 (Weak CD{K,N) spaces have the same growth as then' 
discretizations). Let {X,d,u) be a weak CD{K,N) space, K < 0, N < oo, 
satisfying the following non- collapsing condition: 

(*) There exist tq, Vq > such that V(x,ro) > Vo, for all x ^ X. 

Let G be a discretization of X. Then X has polynomial (exponential) 
volume growth iff G has polynomial (exponential) volume growth. 

As in the case of the proof of Theorem 5.5, the proof closely mimics 
(modulo the necessary modifications) the one of the original result of Kanai, 
as is given in [9], pp. 198-199. 

Proof. First note that, for any r > 0, there exists c(r) > such that 

V(x, r) > c(r) , 

for any x G X. 

Indeed, if r > tq then c(r) = Vo- If r < rg, then by Theorem 13.221 we 
have 



, . f^S^{t)dt , , \lS^{t)dt 

Now, if X G G and y € -B(x, r), by the maximality of the vertex set of G, 
it follows that there exists p ^ M r\ B{y,e), hence d{x,p) < r + e, and it 
follows that 

B{x,r)C U B{p,e). 

peJ^nB(x,r+e) 

Therefore we have 



rr 

V{x,r)< S^it)dt\Afn{x,r + e) 
Jo 



But G is roughly isometric to X, by Theorem 5.5, that is d(pi,p2) < 
ad{pi,p2) + b. Therefore, 

V(x,r) < ci\B{x,ad{pi,p2) + b)\ , 

where B denotes the ball in the combinatorial metric of G and ci is a con- 
stant. 

This concludes the "only if" part of the proof. 
Conversely, for any x (z G and p > 0, we have 

C2\GnB{x,p)\< V{p,e/2) <V{x,e/2 + p), 

p&Gr\B{x,p) 

where C2 = C2(e/2) is a constant. 

But (|5.ip implies that f3{x,p) C B{x,2ep), therefore 
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C2W{x,p)\ < V{x,2e + 2ep) 



which, combined with the previous inequahty concludes the "if" part of the 



Remark 5.12. In [9J, (*) appears as a condition in the statement of the 
Theorem and we followed this approach. For smooth manifolds, this repre- 
sents, however, a consequence of negative Ricci curvature. More precisely, 
VolB{x,r) > Vbr", where Vq = Vo{n), for any r < InjRad(M'")/2. (It ap- 
pears that this is a result due to Croke |12j.) This fact further emphasize 
the need (already underlined in Remark 4.17 and Question 1) for finding 
curvature bounds for the injectivity (conjugacy) radius in metric measure 
spaces. 

Corollary 5.13. Let Xi, X2 be weak CD{K,N) spaces, K < 0,N < 00, 
satisfying condition (*) above. Then, if Xi,X2 are roughly isometric, then 
they have the same volume growth type. 

Proof. Let Gi be a discretization of Xi, i = 1,2. Then, by Theorem 5.4, Gi 
is roughly isometric to Xi, i = 1,2. Since rough isometry is an equivalence 
relation, it follows, that Gi and G2 are roughly isometric, hence they have 
the same growth rate, by Lemma 5.7. Moreover, by Theorem 5.8, Gi has 
the growth rate Xi. □ 



Remark 5.14. By [44j, Theorem 29.9 the results above hold, of course, for 

smooth metric measure spaces as well. 



We have shown that both the triangulation method (and, in consequence, 
the existence of quasimeromorphic mappings) and the discretization results 
extend easily to a wide class of spaces, provided that they satisfy a (gener- 
alized) Gromov-Bishop Theorem. 

Of course, there are also the hidden assumptions that the base space X is 
"rich" enough and the generalized curvature is sufficiently "smooth" . In the 
absence of these conditions, it is not certainly that a Gromov-Bishop type 
theorem can be obtained. Indeed, for the discretization due to J. Ollivier 
[30j . it is not even clear what should be the "reference volume" (and even if 
such a volume exists). Also, unfortunately, this appears to be the case for 
the rough curvature bounds [5]. 
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proof. 



□ 




6. Final Remarks 
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